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This function is originated from the exiled Romanian professor 
Florentin Smarandache . It is defined as follows: 



For any non-null integers n, S (n) is the smallest integer such 
that (S (n) ) ! is divisible by n. 

The importance of the notion is that it characterizes a prime 
number, i . e . : 

Let p > 4, then: p Is prime if and only if S (p) = p. 

Another properties : 

If (a,b) — 1 , then S (ah) = max { S (a) , 5 (b) }; 

and 

For any non-null integers , S (ab) < S (a) + S (b) . 

{All three found and proved by the author in 197 9 {see [3] , 15, 12- 
13, 65) . } 

If n > 1, then S (n) and n have a proper common divisor. 
{Found and proved by student Prodinescu in 1993: as a lemma needed 
to solve the conjecture formulated by the author in 1979 that: 
the equation S (n) = S (n + 1) has no solutions 
(see [3] , 37, and [30] ) . } 

Etc . 



Also, an infinity of open/unsolved problems, involving this 
function, provoked mathematicians around the world to study it and 
its applications (computational mathematics, simulation, quantum 
theory, etc . ) . 

Thus, the unsolved question: 

n 



Calculate lim l 1 + 



- log S (n) J , 



(see [3] , 29) 



n-oo k=2 S (k) 

made by the author in 1979, has been separately proved by J. 
Thompson from USA in 1992 (see [18] , 1) , by Nigel Backhouse from 

United Kingdom in 1993 (see [25] ) , and by Pal Gr0nas from Norway in 
1993 (see [51] ) that this limit is equal to «». 

The author wonderred if it's posible to approach the function 
(see [3] , 197 9, 25-6) , but Ian Parberry expressed that one can 

immediately find an algorithm that computes S (n) in 
0 (nlogn/loglogn) time (see [38], 1993). 

Some unsolved (by now!) other problems stated by the author in 
1979 (see [3] , 27-30) : 

a) To find a general form of the continued fraction expansion 
of S (n) /n, for all n>2 . 

b) What is the smallest k such that for any integer n at 

least one of the numbers S (n) , S (n+1) , ..., S (n+k-1) is a 




perfect square? 

c) To build the largest arithmetical progression a lr a 2 , . . . , 

a r for which their images by the function are also an 
arithmetical progression. 

Etc. 

In 1975 Smarandache was a student at the University of Craiova, 
and he was attracted by the Number Theory. He created and 
published a lot of proposed problems of mathematics in various 
scientific journals. He liked to play with the n umb ers ... 

Thus, in 1980 his research paper "A Function in the Number Theory", 
based on a special representation of integers, was published (for 
the first time) in <Analele Universit&pii Timi?oara>, Seria §tiinpe 
Matematice, Vol. 18, pp . 79-88, 

and was reviewed in <Zentralblatt fur Mathematik>, 471.10004, 1982, 
by P. Kiss, and in the Mathematical Reviews>, 83c: 10008, 1983, by 
R. Meyer. 

In 1988 he escaped from the Ceau§escu' s dictatorship, spent 
almost two years in a political refugee camp in Turkey (Istambul 
and Ankara) , and finally emigrated to the United States . 

Articles, notes, quickies, comments, proposals related to the 
Smarandache Function were presented to international conferences 
within the Mathematical Association of America or the American 
Mathematical Society at the New Mexico State University (Las 
Cruces) , New Mexico Tech. (Socorro) , University of Arizona 
(Tucson) , University of San Antonio, University of Victoria 
(Canada) etc. or published in <Octogon> (Sacele) , <Gazeta 
Matematic3.> (Bucharest) , <The Mathematical Spectrum> (UK) , 
<Elemente der Mathematik> (Switzerland) , <The Fibonacci Quarterly> 
(USA) etc. 

In 1992 Dr. J. R. Sutton from United Kingdom designed a BASIC 
PROCedure to calculate S (n) for all powers of a prime number up to 
a maximum. (see [26] ) 

Jim Duncan from United Kingdom computed up to S(1499999), the 
first million taking 50 hours in Lattice C on an Atari 1040ST. 
(see [17]) 

Also, John McCarthy from United Kingdom estimated that his 
machine would take several years to just calculate and store S (n) 
to disk for the entire range of n it can handle (0<n<2 / '32) , and 
using the compression detailed in ncld9207.c at least 12 Gigabytes 
of disk space would be needed. It took about 3 hours for his 
program to work out that 3,303,302 pages (!) would be needed to 
list the full range of n and S (n) . (see [15]) 

In 1993 Henry Ibstedt from Sweden used a dtk-computer with 
486/33MHz processor in Borland's Turbo Basic and calculated S (n) 
for n upto 10 s which took 2 hours and 50 minutes! (see [52]) 

A group of professors (V. Seleacu, C. Dumitrescu, L. Tupescu, I. 
P&trascu, M. Mocanu) and scientific students from the University of 
Craiova, having a weekly meeting, are doing research on the 
function and its applicability. 
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The Smarandache Function, together with a sample of The Infinity of Unsolved Problems 
associated with it, presented by Mike Mudge . 1 



Smarandache Function, S(n) 
(originated by Florentin Smarandache 
— Smarandache Function Journal , vol 
1. no 1, December 1990. ISSN 1053- 
4792) is defined for all non-null inte- 
gers, n, to be the smallest integer such 
that (S(n))! is divisible by n. 

Note N! denotes the factorial func- 
tion, N!=lx2x3x...xN: for all positive 
integer N. In addition 0! = 1 by defini- 
tion. 

S(n) is an even function. That is, S(n) 
= S(-n) since if (S(n))! is divisible by n 
it is also divisible by -n. 

S(p) = p when p is a prime number, 
since no factorial less than p! has a 
factor p in this case where p is prime. 

The values of S(n) in Fig 1 are easily 
verified. For example, S(14) = 7 be- 
cause 7 is the smallest number such 
that 7! is divisible by 14. 

Problem (i) Design and implement an 
algorithm to generate and store/tabu- 
late S(n) as a function of n. 

Hint It may be advantageous to con- 
sider the STANDARD FORM of n, viz n 

= ep*i p * 2 p*r where e = ± 1 , p, ,p r ...p. 

denote the distinct prime factors of n 
and a., a, a r are their respective mul- 

tiplicities. 

Problem (ii) Investigate those sets of 
consecutive integers (i,i+ l.i+ 2, ... i-*-x) for 
which S generates a monotonic increas- 
ing (or indeed monotonic decreasing) 
sequence. 

Note For (1,2. 3,4 ,5) S generates the 
monotonic increasing sequence 
0, 2,3,4, 5; here i = 1 & x = 4. 

If possible estimate the largest value 
of x. 

Problem (ill) Investigate the existence 
of integers m,n,p,q & k with n*m and 
for which: 

either (A): S(m) + S(m+1) +....S(m+p) = 
S(n) + S(m+1) +....S(n+q) or (B): 

S(mF + S(m+1) 2 +....S(m+p) 2 

S(n) 2 + S(n+i^'' +....S(n+p) 2 = * 

Problem (iv) Find the smallest integer k 



for which it is true that for all n less than 
some given n 0 at least one of: 

S(n), S(n+l)....S(n+k-l) is: 

A) a perfect square 

B) a divisor of k n 

C) a factorial of a positive integer. 

Conjecture what happens to k as n 0 

tends to infinity: i.e. becomes larger 
and larger. 

Problem (v ) Construct prime n umbers 
of the fo rm S (nj S ( n + 1 ) . S (n ■ + k ) : where 
abcdefg denotes the integer formed by 
the concatenation of a.b,c.d ,e,f & g. For 
example, trivially S(2) S(3) = 23 which 
is prime, but no so triv ially 
S(14)S{15)S(16)S(17) = 75617, also 
prime! 

Definition An A- SEQUENCE is an 
integer sequence with l<a i <a 2 <... 

such that no a t is the sum of distinct 
members of the sequence (other than 
aj. 

Problem (vi) Investigate the construc- 
tion of A-SEQUENCES a r a 2 .... such 
that the associated sequences 
S(aJ,S(a 2 ),...are also A-SEQUENCES. 

De finition The k* order forward fi- 
nite differences of the Smarandache 
function axe defined thus: 

D s (xj = *modulus(S(x-i-i) - S(x)), 

Df (x) = D(D(...k-timesD s (x)...)) 
Problem (vii) Investigate the conjec- 
ture that D^ J (l) = 1 or 0 for all k greater 
than or equal to 2. 

c.f. Gilbreath’s conjecture on prime 
numbers, discussed in ‘Numbers Count’ 
PCW Dec 1983. * Here modulus is 
taken to mean the absolute value of 
(ABS.), modulus (y) = y if y is positive 
and modulus (v) = -v if y is negative. 

The following selection of 
Diophantine Equations (i.e. solutions 
are sought in integer values of x) are 
taken from the Smarandache Journal 
and make up: 

Problem (viii) If m & n are given inte- 
gers, solve each of: 

a) S(x) = S(x+1), conjectured to have no 



solution 

b) S(mx+n) = x 

c) Sfmx+n) = m-nx 

d) S(mx+n) = x! 

e) S(x m ) = x n 

fl S(x) m = S(x n ) 

g) S(x) + v = x + S(v), x&v not prime 

h) S(x) + S(v) = S (x+y) 

i) S(x-ry) = S(x)S(y) 

j) S(xv) = S(x)S(y) 



Review, July 1992 

The Smarandache Function: 

a first visit? 2 

This topic is certain to be revisited in 
the near future, and the lack of space 
available here will certainly be rem- 
edied on that occasion. Suffice it to 
report that Jim Duncan computed up to 
S(l499999). the first million taking 50 
hours in Lattice C on an Atari 1040ST. 
In Problem (ii), no evidence for a largest 
value ofx was found, while in Problem 
(vii) the conjecture was verified for the 
first 32,000 values of S(n). The very 
worthy prizewinner is John McCarthy 
of 17 Mount Street, Mansfield, Notts 
NG19 7 AT, who has extensively inves- 
tigated the computation of S(n) up to 
2 32 ; arriving at conclusions such as: 
‘several years of computing', ‘at least 
12Gb of disk space’ and ‘3,303,302 pages 
of output’. John’s concluding comment, 
‘Am I mad?’, is clearlv answered NO! 
by examining his specimen pages of 
output including those relating to 10- 
digit values of n. Listings supplied. 
Details from John directly upon request. 




1 

Republished from <Personal Computer fforld>, No. 112, 420, July 
1992 (with the author permission) , because some of the following 
research papers are referring to these open problems. 

2 

Republished from <Personal Computer World>, No. 117, 412, December 
1992 (with the author permission) . 




ALGORITHM IN LATTICE C TO GENERATE S (n) 



by Jim Duncan 

9 Ryeground Lane 

Formby 

Liverpool 

L37 7EG ENGLAND 



Computer: Atari 1040ST 

Run time: ca, 50 hrs to generate S(n) for 1 000 000 numbers 



#i nclude <stdio.h> 

^include <math.h> 

unsigned long int pst,ast,s,t; 



ma i n < ) 
i 

long int n; 

FILE *fpy 

print* ("input n (1 <- n < 2 147 483 647)\n">; 
scant ("Xld" , tin); 

■fp * * open < "PRN: M , "w" ); 

tpr intt(*p, "n * XI d S(n) =* 7.1 d\n" , n, smaran(n))| 
f cl ose (f p ) ; 

> 

smaran <m) 

unsigned long int m; 

< 

unsigned long int <nst , p , a , * act ; 
double r ; 

H ( ro == i) 

return (0) ; 

else < /* STANDARD FORM ot m */ 

r = mst 33 m ; 
p = i; 
tact =* 1; 

while < ++p <= sqrt(r)) i 
a =» 0? 

while (mst X p == 0) { 

mst - mst/p; 
a++ ; 

> 

r =* mst; 
it (a > 0) < 

pst = p; 
ast « a; 
t * s » 0; 
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(t) with */ 



tors3 ( ) ; /* 'find smallest factorial 
/* p^a divsor */ 
if (t > fact) 

fact * t; 

> 

> 

if (mst > fact) 

fact = mst; 
return (fact) ; 

> 

> 

torsi < ) 

C 

unsigned long i nt i; 

/* test number is pst^ast */ 

/* s is the difference between a factorial number t and ast* {pst- 1 ) */ 

/* s forms a pattern which determines the smallest value of t for which */ 
/* the test number is a divisor */ 

i = 0; 

while ( ++i < pst*pst && t-s < ast*(pst-l>) { 
if <i X pst == 0) 

s = s-pst+2; 

el se 

s++; 

t += pst; 

> 

> 



tor s2 ( ) 

< 

unsigned long int i; 

torsi < > ; 
i = 0; 

while ( ++i < pst*pst && t-s < ast* <pst-l ) ) < 

if (i X pst == 0) 

s = s-3*pst+4; 

el se 

s = s-2*pst+3; 
t +* pst; 
torsi < ) ; 



tors3 ( ) 

< 

unsigned long int i; 

tors2 ( ) ; 
i = 0; 

while (++i < pst*pst && t-s < ast*(pst-l)) { 

if (i X pst == 0) 

s = s-5*pst+6; 

el se 

s = s-4*pst+5; 
t +=■ pst; 
tors2 ( ) ; 

> 
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MONOTONIC INCREASING AND DECREASING SEQUENCES OF S (n) 

by Jim Duncan 



Prob 1 efn ( 11 ) 



Monotonic increasing and monotonic decreasing sequences of S(n) 
were investigated for x >= 6. 



First number 


(i ) 






number of 


sequences S 






in range 






X 


= 6 


x = 


7 


x = 8 


X 


= 9 








i nc 


dec 


inc 


dec 


• i nc dec 


i nc 


dec 


1 


499 


999 


75 


83 


7 


10 


0 2 


0 


0 


500 000 


999 


999 


80 


76 


14 


18 


1 3 


1 


1 


1 000 000 - 1 


499 


999 


75 


63 


8 


10 


1 2 


1 


1 



There appears to be no evidence for a largest value for x. The sequences 
for x = 9 are shown in Results Table 1. 



The existence of sequences with the same first order finite differences 
was then considered eg: 



440 


S ( i ) 


= 11; 


i = 5073 


S ( i ) 


89 


441 


S ( i +1 ) 


= 14; 


i + 1 = 5074 


S ( i + 1 ) * 


59 


442 


S(i+2) 


= 17 ; 


i+2 = 5075 


S(i+2> = 


29 



Apart from the initial quartet 2,3,4, 5 all such sequences with 
i K 1 000 000 are triplets. If the first order finite differences 
are multiples of 6 then the S(n) values appear to be prime numbers. 
The values are shown in Results Table 2. 
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RESULTS 



TABLE 1 



Sequences of S (n) 



X 


— 


9 








n 


= 


586951 


S < n ) 


=: 


586951 


n 


= 


586952 


S (n) 


= 


73369 


n 


= 


586953 


S (n) 


= 


21739 


n 


= 


586954 


S (n) 


= 


9467 


n 


- 


586955 


S (n) 


= 


1319 


n 




586956 


S ( n ) 




1 193 


n 


= 


586957 


S(n) 


= 


1181 


n 


= 


586958 


S < n ) 


= 


1091 


n 


= 


586959 


S (n) 


= 


677 


n 




586960 


S < n ) 


= 


29 



X 


S3 


9 








n 


23 


721970 


S (n) 


= 


73 


n 


= 


721971 


S (n) 


= 


827 


n 


= 


721972 


S <n) 


= 


907 


n 




721973 


S < n ) 


=3 


6067 


n 


= 


721974 


S <n ) 


= 


10939 


n 




721975 


S(n> 


=3 


28879 


n 


= 


721976 


S <n ) 


as 


90247 


n 


CS 


721977 


S <n ) 


= 


240659 


n 


= 


721978 


S (n) 


= 


360989 


n 




721979 


S (n ) 




721979 



x =* 9 



n 


= 


1091 150 


S <n ) 


= 


157 


n 


- 


1091151 


S (n) 




709 


n 


= 


1091152 


S ( n ) 


= 


1451 


n 


= 


1091153 


S ( n ) 


23 


1607 


n 


2= 


1091154 


S <n ) 


= 


6271 


n 


3 


1091155 


S (n) 




16787 


n 




1091156 


S(n) 




24799 


n 


= 


1091157 


S(n) 


= 


363719 


n 


= 


109 J 158 


S ( n ) 


as 


545579 


n 




1091159 


S (n) 


= 


1091159 



X = 9 



n 




1473257 


S (n ) 


= 


1473257 


n 


=3 


1473258 


S < n ) 


S3 


8467 


n 


= 


1473259 


S ( n ) 


= 


6323 


n 


3 


1473260 


S (n) 


3 


3877 


n 


= 


1473261 


S Cn) 


= 


3533 


n 


3 


1473262 


S(n) 


=3 


2239 


n 


3 


1473263 


S(n) 


= 


1999 


n 


3 


1473264 


S ( n ) 


=* 


787 


n 


= 


1473265 


S(n) 


33 


557 


n 


3 


1473266 


S ( n ) 


- 


463 




RESULTS 



TABLE 2 



same difference = 1 n = 4 S(n- 2 ) = 2 S(n-l) = 3 S(n) = 4 

same difference =1 n = 5 S(n-2) = 3 S<n- 1 ) = 4 S<n) = 5 

same difference =11 n = IB S(n-2) = 6 S(n-l) = 17 Sin) = 6 

same difference =3 n = 442 Sin-2) = 11 S(n-l) = 14 S(n) = 17 

same difference =30 n = 5075 Sin-2) = 89 S(n-l) = 59 Sin) = 29 

same difference = 60 n = 6409 Sin-2) = 149 S(n-l) = 8 Q S(n> = 29 

same difference =40 n = 6478 S(n-2> = 127 S(n-l) = 79 3(n) = 31 

same difference = 36 n = 8177 Sv'n-2) = 109 S(n-l) = 73 S(n> = 37 

same difference =84 n = 13717 S(n-2! = 211 S(n-l) = 127 Sin' - 43 

same difference = 168 n = 20468 S(n-2) = 379 S(n-l) = 211 Sin) = 43 

same difference = 210 n = 22591 S(n-2) = 461 Sin- 1 ) = 251 Sin) = 41 

same difference = 120 n = 35145 S(n-2) = 311 Sin- 1 ) = 191 Sin) = 71 

same difference = 180 n = 59719 Sin-2) = 449 Stn- 1 ) = 269 Sin) = 89 

same difference = 150 n = 67771 Sin-2) = 401 Sin- 1 ) = 251 Sin) = 101 

same difference = 264 n = 73425 Sin-2) = 617 S(n- 1 > = 353 Sin) = 89 

same difference = 840 n = 74005 Sin-2) = 1721 Sin- 1 ) = 881 Sin) = 41 

same difference =24 n = 82297 Sin-2) = 151 Sin- 1 ) = 127 Sin) = 103 

same difference =60 n = 104669 Sin-2) = 251 Sin- 1 ) = 191 Sin) = 131 
same difference = 330 n = 111507 Sin-2) = 769 Sin- 1 ) = 439 Sin) = 109 
same difference =36 n = 114427 Sin-2) = 199 Sin- 1 ) = 163 Sin) = 127 
same difference = 252 n = 120523 S(n-2) = 631 Sin- 1 ) = 379 Sin) = 127 

same difference = 120 n = 129928 Sin-2) = 389 Sin- 1 ) = 269 Sin) = 149 

same difference = 952 n = 146004 Sin-2) = 1973 Sin- 1 ) = 1021 Sin) = 69 

same difference = 600 n = 153520 Sin-2) = 1301 Sin- 1 ) = 701 Sin) = 101 



same 


di i i erence 




12 


n 


= 180482 S (n-2) 


= 47 S(n-l) = 


59 Sin) = 71 


same 


di -f-ference 


=S 


660 


n 


= 181485 S (n-2) 


= 1429 S(n-l) 


* 769 S < n ) 


=» 109 


same 


di f-ference 


= 


60 


n 


- 189954 S (n-2) 


= 53 S (n— 1 ) = 


113 S ( n ) = 


173 


same 


di fference 


as 


90 


n 


* 192067 S (n-2) 


= 359 S(n-l) = 


= 269 S ( n ) = 


= 179 


same 


d i f-ference 


= 


324 


n 


= 198697 S (n-2) 


* 811 S(n-l) 


= 487 S ( n > 


= 163 


same 


di Herence 


= 


336 


n 


= 209752 S (n-2) 


= 839 S(n-l) 


= 503 S (n ) 


= 167 


same 


di i i erence 


= 


228 


n 


= 227099 S (n— 2) 


» 647 S(n-l) 


= 419 S (n ) 


= 191 


same 


d i Herence 


* 


150 


n 


= 231039 S (n-2) 


= 499 S (n— 1 ) 


= 349 S ( n ) 


= 199 


same 


di -f *f erence 




264 


n 


= 253725 S (n-2) 


= 727 S(n-l) 


= 463 S(n) 


= 199 


same 


di f i erence 


= 


210 


n 


= 266915 S (n-2) 


= 631 S (n— 1 ) 


= 421 S (n ) 


= 21 1 


same 


di i erence 


= 


648 


n 


= 297638 S (n-2) 


= 1459 S(n-l) 


= 811 Sin) 


= 163 


same 


di i i erence 


= 


2808 


n = 306128 Sin-2 


) = 5669 Sin-1) = 2861 Sin) = 53 


same 


di -f *f erence 


= 


1320 


1 


n = 324384 Sin-2) = 2749 Sin-1) = 1429 Sin) = 109 


same 


di ii erence 


■ 


18 


n j 


= 326163 Sin-2) 


= 199 Sin-1) = 


131 S (n ) = 


163 


same 


di Terence 


= 


240 


n 


= 342965 Sin-2) 


= 719 Sin-1) 


- 479 S ( n ) 


= 239 


same 


di i i er ence 


— 


36 


n : 


= 346390 Sin-2) 


= 139 Sin-1) = 


103 S < n ) = 


67 


same 


di f f erence 


= 


300 


n 


= 386906 Sin-2) 


= 47 Sin-1) = 


347 5 (n ) = 


647 


same 


di Terence 




840 


n 


= 409422 Sin-2) 


= 1861 Sin-1) 


= 1021 Sin) = 181 


same 


di -f -f erence 


= 


270 


n 


= 440375 Sin-2) 


= 811 Sin-1) 


= 541 S (n ) 


= 271 


same 


d i Herence 




936 


n 


= 443450 Sin-2) 


= 2053 Sin-1) 


= 1117 Sin) = 181 


same 


di f -F erence 


* 


120 


n 


= 443850 Sin-2) 


= 509 Sin-1) 


= 389 S ( n ) 


= 269 


same 


d i -ff erence 


= 


792 


n 


= 443969 Sin-2) 


= 1783 S in- 1 ) 


= 991 S ( n ) 


= 199 


same 


di -f-ference 


= 


450 


n 


= 450043 Sin-2) 


= 1151 Sin-1) 


= 701 S <n ) 


= 251 


same 


d i ii erence 


= 


306 


n 


= 451215 Sin-2) 


= 883 Sin-1) 


= 577 S < n ) 


= 271 


same 


di f -f erence 


= 


210 


n 


= 460559 Sin-2) 


= 701 Sin-1) < 


= 491 S ( n ) 


= 281 


same 


di f ference 


= 


240 


n 


= 464212 Sin-2) 


= 761 Sin-1) 


= 521 S (n ) 


= 281 


same 


di fference 


= 


360 


n 


= 470727 S in— 2) 


= 991 Sin-1) = 


= 631 S (n ) 


= 271 


same 


di i i erence 


= 


624 


n 


= 473922 Sin-2) 


= 1481 Sin-1) 


= 857 S ( n ) 


= 233 


same 


di i 4 erence 




90 n = 


= 481779 Sin-2) = 


= 449 Sin-1) = 


359 S ( n ) = 


269 


same 


di i i erence 


= 


126 


n 


= 511688 Sin-2) 


= 131 Sin-1) = 


= 257 S(n) = 


= 383 


same 


d i -f f er ence 


= 


672 


n 


= 512894 Sin-2) 


= 1583 Sin-1) 


= 911 S (n) 


= 239 


same 


di -f i erence 


= 


480 


n 


= 521946 Sin-2) 


= 1231 Sin-1) 


= 751 S ( n ) 


= 271 


same 


d i *f i erence 


= 


714 


n 


= 531775 S (n-2) 


= 1667 Sin-1) 


= 953 S (n ) 


= 239 


same 


d i -f *f erence 


- 


726 


n 


= 543455 Sin-2) 


= 1693 Sin-1) 


= 967 8 ( n ) 


= 241 




same difference = 306 n = 565187 S(n-2) = 919 S(n-l) = 613 S(n) = 307 

same difference = 552 n » 574498 S(n-2> - 1381 S(n-l) = 829 S(n) = 277 

same difference = 1176 n = 586272 Sin-2) = 2549 S(n-l) = 1373 S(n) » 197 

same difference = 3444 n = 592537 S(n-2) = 6971 S(n-l) = 3527 S(n> = 83 

same difference = 390 n = 609871 S(n-2> = 1091 S(n-l) = 701 S(n) = 311 

same difference = 840 n = 629508 S<n-2) = 1931 S(n-l) = 1091 S(n) = 251 

same difference = 336 n = 681077 S(n-2) = 1009 Sin-1) = 673 S(n) = 337 

same difference = 612 n = 705793 S(n-2) = 1531 S(n-l) = 919 S(n) = 307 

same difference = 78 n = 724319 S(n-2) = 467 S(n-l) = 389 S(n) = 311 

same difference = 264 n = 726827 S(n-2) = 881 S(n-l) = 617 S(n) = 353 

same difference = 498 n = 731179 S(n-2> = 1327 S(n-1> = 829 Sin) = 331 

same difference - 240 n = 751746 Sin-2) = 839 Sin-1) = 599 Sin) = 359 

same difference = 1356 n = 778837 S<n-2> = 2939 S(n-l) = 1583 Sin) '=227 
same difference = 1020 n = 792675 S(n-2) = 2311 Sin-1) = 1291 Sin) = 271 

same difference = 2214 n = 803427 Sin-2) = 4591 Sin-1) = 2377 Sin) = 163 

same difference = 1590 n = 810451 Sin-2) = 3391 Sin-1) = 1801 Sin) = 211 

same difference = 252 n = 837969 Sin-2) = 883 Sin-1) = 631 Sin) = 379~ 

same difference = 552 n = 898783 Sin-2) = 1471 Sin-1) = 919 Sin) = 367 

same difference = 2850 n = 930311 Sin-2) = 5851 Sin-1) = 3001 Sin) = 151 
same difference = 540 n = 941057 S(n-2) = 1459 Sin-1) = 919 Sin) = 379 
same difference = 240 n = 943553 Sin-2) = 881 Sin-1) = 641 Sin) = 401 

same difference = 120 n = 975546 Sin-2) = 619 Sin-1) = 499 Sin) = 379 

same difference = 1122 n = 997443 Sin-2) = 2551 Stn-1) = 1429 Sin) = 307 
same difference = 264 n = 1026550 Sin-2) = 947 Sin-1) = 683 Sin) = 419 
same difference = 684 n = 1028985 Sin-2) = 1747 Sin-1) = 1063 Sin) = 379 

same difference = 744 n = 1042162 Sin-2) = 1861 Sin-1) = 1117 Sin) = 373 

same difference = 510 n = 1053175 Sin-2) = 1429 Sin-1) = 919 Sin) = 409 
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ON THE CONJECTURE D, w (l)=l or 0 for k>=2 

by Jim Duncan 



Pr obi em (vi 1 ) 

For the first 32 000 S(n) 's the conjecture that D s <*> (i) = 1 or 0 for 

k 2 true * The ratio of the number of ones to the number of zeros 
appears to be approximately 1 for large values of k. The results are 
shown in Results Table 3. 



The true differences S(x-l) - S(x) 
differences D S *“«><1> were found t 
For large values of k ( > 100) the 
approximately equal to -2. Some va 



RESULTS TABLE 3 

Ratio: 

k D 3 (k > (1) total l's 

total 0 ' s 



were calculated and the k th order 
i increase rapidly with increasing k. 
ratio D, •“*> (1) (l) i s 
ues are shown in Results Table 4. 



D 1 =2 
D 1 000 = 


0 


stl/sto 


— 


1.11 205 1 


D2000 = 


0 


stl/stO 


= 


1 . 056584 


D3000 = 


0 


st 1 /stO 


= 


1 . 075433 


B4000 = 


0 


stl/stO 


= 


1 . 015625 


D50O0 = 


0 


st 1 /stO 




1.01 0S6 1 


D6000 = 


1 


st 1 /stO 


= 


0. 991039 


D7000 = 


0 


stl/stO 


- 


0. 990048 


DB000 = 


0 


st 1 /stO 


= 


0. 977014 


D9000 = 


0 


stl/stO 


= 


0. 987412 


D 10000 = 


0 


stl/stO 


= 


0. 998201 


D11000 = 


1 


stl/stO 


= 


1.011154 


D 12000 = 


1 


stl/stO 


= 


1 . 010556 


D 13000 = 


1 


stl/stO 


= 


1 . 015036 


D 14000 = 


0 


stl/stO 


= 


1 . 018601 


D 15000 = 


1 


stl/stO 


= 


1.012748 


D 16000 = 


0 


stl/stO 




1 . 004134 


D 17000 = 


0 


stl/stO 


= 


1 . 005308 


D 18000 = 


0 


stl/stO 


= 


1 . 004789 


D 19000 = 


1 


stl/stO 


= 


1 . 003903 


D20000 = 


1 


stl/stO 


= 


1 . 00471 1 


D21000 = 


0 


stl/stO 


= 


1 . 008129 


D22000 = 


0 


st 1 /stO 


— 


1 . 004830 


D23000 = 


0 


stl/stO 


= 


1.004620 


D24000 « 


0 


stl/stO 


= 


1 . 003590 


D25000 = 


1 


stl/stO 


= 


1 . 004571 


D26000 = 


1 


stl/stO 


= 


1 . 001001 


D27000 = 


0 


stl/stO 


= 


1 . 001260 


D28000 = 


1 


st 1 /stO 


= 


1 . 004080 


D29000 = 


0 


stl/stO 


= 


1 . 006018 


D30000 = 


1 


st 1 /stO 


= 


1 . 0054 1 5 


D31000 = 


0 


stl/stO 


= 


1 . 006408 


D32000 = 


1 


st 1 /stO 


= 


1 . 004699 




RESULTS 



TABLE 4 





k D/ (IO (l) 




1 (i) 






d 3 


* (k-i) ( i ) 


0951 


= 2 . 24442 1 E+23S 


rat i o 


= -2 - 00297 4 


0952 


= -4. 496445E+288 


rat i o 


= -2.003387 


0953 


= 9 . 0099 1 6E+288 


ratio 


= -2.003787 


D954 


= — 1 . 8057 40E+289 


r at i o 


= -2.004170 


0955 


= 3.619671 E+289 


ratio 


= -2.004535 


D956 


= -7. 257009E+289 


r at i a 


= -2.004881 


0957 


= 1 . 455178E+290 


rat i o 


= -2.005204 


0958 


= —2. 918366E+290 


rat i o 


= -2.005504 


0939 


= 5. 853595E+290 


rati o 


= -2.005778 


0960 


= -1 . 174247E+2^1 


rat i o 


= -2.006026 


0961 


= 2. 355827E+291 


ratio 


= -2.006246 


0962 


= -4. 7268 1 9E-t-29 1 


rati o 


= -2.006437 


0963 


= 9. 484823E+291 


ratio 


= -2.006588 


D964 


= -1 . 903346E+292 


rat i o 


= —2 . 006723 


0963 


= 3.81 9685E+292 


ratio 


= -2.006Q27 


0°66 


= —7. 665702E+292 


rat i o 


= —2 . 006894 


0967 


= 1 . 538452E+293 


rat i o 


= -2.006929 


0968 


= —3. 087569E+293 


rati o 


= -2.006932 


D969 


= 6. 196454E+293 


r a t i o 


= -2.006904 


D970 


= -1 . 243531E+294 


rat i o 


= -2.006843 


0971 


= 2. 495458E+2 <3 4 


ratio 


= -2.006751 


0972 


= —5. 007457E+294 


rat i o 


= -2.006629 


D973 


= 1 . 004734E+295 


ratio 


= -2.006476 


0974 


= -2.015791E+295 


rati o 


= -2.006293 


D975 


= 4. 043841E+295 


rati o 


= -2.006081 


D976 


= -8. 111306E+295 


rati o 


= -2.005842 


0977 


= 1.626784E+296 


rati o 


= -2.005576 


D97S 


= —3. 262162E+296 


rati o 


= -2.005283 


0979 


= 6. 540525E+296 


rati o 


= -2.004966 


D980 


= -1 .31 1 130E+297 


rati o 


= -2.004625 


0981 


= 2. 627848E+297 


rati o 


= -2.004262 


0932 


= -5. 265885E+297 


ratio 


= -2.003877 


0983 


= 1.055006E+298 


ratio 


= -2.003473 


D984 


= -2. 1 13228E+298 


ratio 


* -2.003049 


D985 


= 4. 23 1 969E+298 


rati o 


= -2.002608 


D986 


= —8. 473041E+298 


ratio 


= -2.002151 


0987 


= 1 . 69603 1E+299 


ratio 


= -2.001679 


0988 


= —3. 394086E+299 


ratio 


= -2.001 193 


0989 


= 6. 790531E+299 


ratio 


= -2.000695 


0990 


= -1 . 358233E+300 


rat i o 


= -2.000186 


0991 


= 2. 716013E+300 


rati o 


= -1.999667 


0992 


= -5. 429688E+300 


ra+* i o 


= -1.999139 


D993 


= 1 . 035 1 SOE+ 30 1 


ratio 


= - 1 . 998604 


0994 


= -2. 16S257E+301 


ratio 


= -1.998063 


0995 


= 4.331 129E+301 


rat i o 


= -1.997516 


D996 


= -8. 6491 19E+301 


ratio 


= -1.996966 


D997 


= 1 . 726721E+302 


r at i o 


= -1.996413 


0998 


= -3. 44629 1 E+302 


ratio 


= -1.995858 




A Simple Algorithm to Calculate S(n) 

by John C. McCarthy 



Introduction 



This short paper first outlines 
calculating S(n) (the smallest integer m 
n). Doubtless, there exist more subtle 
hope some readers will devise these and 
through this journal. 



an "obvious" algorithm for 
such that m! is divisible by 
and efficient algorithms. I 
enlighten me concerning them 



This is followed by a small scale investigation of the efficiency 
of the algorithm. 

Then there is a short discussion of a simple way of reducing the 
space required for storage of all S(n) for ranges of n. The storage 
space required for S(n) for all n which my routines can handle is 
considered . 

Heavily commented listings of an implementation of the algorithm 
in "C", sample output and timing data are included to help illustrate 
the algorithm. 



The Algorithm 

The algorithm is described in detail at the start of the header 
file "S(n).H". Together with "S(n).C", this forms all the code 
necessary to implement the algorithm. Note that, for the S(n) function 
to work correctly, the function make_pr imes ( ) must first be called from 
the main program. 



The code for printing S(n) and timing the routines has been 
omitted. These activities are both implementation specific and easily 
done. They are therefore left as an exercise for the interested 
reader . 

The algorithm hinges on finding the prime factors of n. 
Improvements on how this is done will most benefit its efficiency. 



To be practical, the given implementation of the algorithm only 
works for 0<n<2 32 . However, the algorithm is generally applicable to 
any non-null integer. 

Tables of S(n), constructed using the routines of "S(n).C", for 
the largest 2000 permitted n are included. 1 My paging 
routines are rather elaborate. Using them (without printing!), it took 
2.4 hours to discover that 3,745,708 pages, as tightly packed as those 
shown, would be required to print S(n) for all 0<n<2 32 . 



Efficiency of the Algorithm 

In a letter to R. Muller (about computing the Smarandache 
Function, July 19, 1993), Ian Parberry (editor of <SIGACT News>, 



2 

For the smallest 4800 numbers, see Ibstedt' a table (pp. 43-50) of 
this current journal. 
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Denton, Texas) expressed that one can immediately find an algorithm 
that computes S(n) in O(nlogn/loglogn) time ('A Brief History of the 
"Smarandache Function"' by Dr. Constantin Dumitrescu, Department of 
Mathematics, University of Craiova, Romania). Disappointingly, a 
little analysis of the accompanying timing data on my TI85 advanced 
scientific calculator reveals that my algorithm is somewhat worse than 
this. 



Trying to fit the version 2 timing data to various 0(f(n)), I 
obtained the following results (x=3355443200 and 10 (0( x+99 ) -0( x- 100 ) ) 
is calculated for comparison with the last entry of the version 1 
t iming data ) .• 



O ( f ( n ) ) 


Correlation 

Coefficient 


0( 2 3 2 -l ) 
( years ) 


10 (0 ( x+99 ) -O ( x- 
(mi 1 1 i seconds 


0(n) 


0.9928879 


0.6092 


8909 


O(nlogn/loglogn) 


0.9944006 


0.7906 


11827 


0(n/n) 


0.9997756 


24.2 


469178 



0(n/n) fits the version 2 timing data best, although the time it 
predicts for the last entry of the version 1 timing data is almost 3 
times too large. Hence, I assume the time complexity of my algorithm 
is a little better than 0(n/n). 

As a rough upper limit on the time my program (on my 20MHz 368DX 
PC) would take to calculate S(n) for all 0<n<2 3Z , let us assume that 
every value of n requires as much time as each n in the range of the 
last entry of the version 1 timing data (= 159111/199/10 » 79.9553 ms). 
In this "worst case", it would take 10.882 years. 0(n/n) time 
complexity predicts more than twice this value, which is a measure of 
how pessimistic it is. 

I would welcome a more rigorous analysis of the time complexity of 
my algorithm as I presently lack the necessary expertise. 



Simple Compression of Stored S(n) 

Without compression, each S(n) would be stored as a 32-bit (= 4 
bytes) value. Hence 2 34 bytes (= 16 Gigabytes) would be required to 

store S(n) for all 0<n<2 3Z . 

This requirement can be reduced considerably if we use the high 
bit of each each byte of each value to indicate if it is the last byte 
of the value. If the bit is set it means that further byte(s) are 
required and if it is reset it means that the byte is the last byte of 
the current value. This means that only 7 bits of each byte are used 
to form the numerical part of the value. Assuming that, as with Intel 
format, the values are stored low-'byte' (actually 7 bits) first, here 
are some examples: 

i) 127 requires seven bits and so just one byte (with high bit 
reset to indicate no further bytes). 

ii) 16,000 requires 14 bits. So it is stored as two bytes. The 
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first is 0 (16,000 mod 128) + 128 (to set the high bit 
indicating there is more to come). The second is 125 (16,000 
div 128) (with high bit reset to indicate no further bytes). 
This reads simply as 0 (with more to follow) + 128*125 (no 
more to follow). 

iii) A number stored as the three bytes 57+128, 93+128 and 125+0 
would similarly represent: 

57 + 93*128 + 125*128*128 = 2,059,961. 



The largest numbers that can be 
bytes is thus as follows: 

1 byte can code up to 2 7 -l = 

2 bytes can code up to 2 14 -1 = 

3 bytes can code up to 2 Z1 -1 = 

4 bytes can code up to 2 Z *-1 = 

5 bytes can code up to 2 3 5 -l = 

largest unsigned long). 



represented by a given number of 



127. 

16,383. 

2,097,151. 

268,435,455. 

34,359,738,355 (or 8 times the 



For small values of n, the savings are considerable (400%). 
However, even large n often have small S(n). 



Using this 


technique to 


compress all 


S(n) 


calculated for some 


ranges of n (each range was 
resul t s : 


also stored) 


, I obtained the following 


range of n 


compression 


time taken 
( seconds ) 


size 


size after pkzip 


1 


wi thout 


4.5 


40,008 


19,836 


-10,000 


wi th 


4.7 


15,749 


15,267 


2,147,478,648 


wi thout 


827.3 


40,008 


33,729 


-2,147,488,647 


wi th 


842.4 


33,541 


30,836 


4,294,957,296 


wi thout 


1,066.2 


40,008 


34 , 320 


-4,294,967,295 


wi th 


1,085.1 


34,330 


31,634 



The results indicate that this compression is a little better than 
pkzip’s (a commercial file compression utility). Application of pkzip 
to a pre-compressed file also gives a slight improvement. 



Assuming that the savings shown for the middle range of 10,000 n 
are the average of all ranges of 10,000 n, using my compression 
together with that of pkzip would permit storage of S(n) for all 
0<n<2 3Z in about 3.0836*2 3Z = 12.3344 Gigabytes. So look out for sets 
of 19 CD-ROMs with all your favourite numbers on them! 



21st November 1993 
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S ( n ) . H 



/* (c). 1993. 11. 13 .John. C. McCarthy 
”S(n) .h" 

Example Implementation of A Simple Algorithm to Calculate S(n), 

The Smarandache Function: 

Because there are more people familiar with C than with C++, this module 
has been written entirely in C (apart from "//" style comments). The 
module was compiled using Borland C++ version 3.1. 

For efficiency, n is constrained to the limits of an unsigned long. Hence, 
0 <= n <= 2~32 - 1 ( = 4,294,967,295). represents exponentiation). 

Although catering for n of vast magnitude is possible, it imposes heavy 
storage and processing overheads. The range of an unsigned long therefore 
seems a reasonable compromise. 

The algorithm depends on the most elementary properties of S(n): 

1) Calculate the STANDARD FORM (SF) of n: 

In SF: n = +/- (pl“al )* (p2“a2 )*...* (pr~ar ) where pi, p2,...pr denote the 
distinct prime factors of n and al, a2,...ar are their respective 
multiplicities. 

2) S(n) = max [ S (pl~al ),..., S(pr~ar )] . 

3) S(p~a), where p is prime, is given by: 

3.1) a<=p ==> S(p~a) = p*a . 

3.2) a>p ==> S(p~a) = x < p*a. In this case, fortunately rare, x 

is the smallest integer such that p appears as a factor in the 
list of all integers > 1 and <= x at least a times. Let the no. 
of times p appears as a factor in the list of all integers > 1 
and <= y be f(y, p). Then: 

f(y, p) = Z tint (y/(p~i ) ) ] for i > 0 while y>=(p“i). 

Hence, x is the smallest integer such that f(x, p)>=a. 

Note that between succesive integer multiples of p there are no 
integers which have p as a factor. The trick here is to look for 
the largest multiple of p (call it c), such that f(p*c, p)<=a 
(so that x = p*c, if f(p*c, p) = a, else x = p*(c+l)): 

3.2.1) c = a-2 (largest possibilty for c since f(p*(a-l), p)>=a 
when a>p (Note: f(p*(a-l), p)=a is not sought for slight 
performance gain)). 

3.2.2) z = f (p*c , p) . 

3.2.3) While(z>a): 

3. 2. 3.1) d = no. of times p appears as a factor of p*c 

= (no. of times p appears as a factor of c) + 1. 

3. 2. 3. 2) c = c-1 (next largest possibility for c). 

3. 2. 3. 3) 2 = z-d (= f(p*c, p) ) . 

3.2.4) If ( 2 <a) , x = p* ( c+ 1 ) . 

3.2.5) Else x= p*c . 

To calculate the prime factors of all 32-bit n requires use only of primes 
< (2~16) (i.e. all primes expressible as an unsigned short integer). This 

is because any factor of n remaining after division of n by all its prime 
factors < (2“16) is simply a prime. Since there are only 6542 16-bit 
primes, the program first creates a list of these (which only takes about 
4 seconds on my 20 MHz 386DX PC) so that they never have to be 
recalculated, thus saving much time. 

*/ 
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S (n) .H 



^define PRIMES16 6542 // The number of 16-bit primes 

#define MAX_SFK 9 /* max. distinct primes in the SF of n. The smallest 

number with more than 9 distinct primes is the product of the 10 smallest 
primes ( = 6,469,693,230), which is substantially more than the largest 
integer expressible as an unsigned long. Hence, 9 distinct primes are 
more than ample. 

*/ 



typedef unsigned long u_long; 
typedef unsigned int u_int; 
typedef enum {false, true} boolean; 

struct SF_struct { 
int s f k ; 

u_long sfp[MAX_SFK] ; 
int s f a [MAX_SFK] ; 

} ; 

extern u_int prime [PRIMES16+1 ] ; // list of all 16-bit primes 

// plus terminating zero. 

void make_primes (void) ; // construct list of all 16-bit primes (primet]). 

// Must be called before calls to getSFU or SO. 
void getSF(u_long n, struct SF_struct *SF); // calc. SF of n and store in SF 
u_long S(u_long n) ; // calc. S(n) 

u_long Spa(u_long p, int a); // calc. S(p“a) where p is prime 
int f(int x, int p) ; I* the number of times the prime p appears as a factor 
in the integers from 1 to x inclusive. This function is only called from 
Spa(p, a) when a>p with x=p*(a-2) (refer to item (3) of algorithm outline 
above). Max value of (a) occurs when p is a minimum, n is a maximum and 
(p~a)=n. So, ( 2“max ( a ) ) =max ( n ) = ( 2“ 32 ) -1 - Hence max(a)<32. So, x<60 
when (a) is at its max. Max value of p (and x) occurs when a=p+l and 
( p~a ) =max ( n ) . So, max (p ) “ (max (p ) + 1 ) = ( 2~ 32 ) - 1 . The upshot is that 
max(p)=9 when a=10. Hence, max(x)=72. This explains why it is safe for 
x, p and the return value of f(x,p) to be passed as ints. 

*/ 



// no. of distinct primes 
// the distinct primes 
// respective multiplicities 



23 




S (n) .C 



/* (c). 1993. 11. 13. John . C . McCarthy 
" S { n ) .c" 



Example Implementation of A Simple Algorithm to Calculate S(n), 

The Smarandache Function: 

This is the code for the module. Refer to ”S(n).h" for details. 

*/ 

#include "S(n) .h” 

u_int prime [PRIMES16+1 ] ; // allocate storage for list of all 16-bit primes 

// plus terminating zero. 

void make_primes ( void) 

{ 

u_int *pp ; // ptr to last prime so far of prime list 

u_int * tp ; // ptr to current test prime 

u_int p; // number being tested for primal ity 

pp=prime; // point to start of prime list 
*pp=2 ; // set first prime to 2 

*++pp=3; // set second prime to 3 

P =5 <‘ // next possible prime. N.B. p is kept odd so that trial 

// division by 2 is unnecessary, 
whi l*e ( true ) { // infinite loop!: 

tp=prime+l; // point to first odd test prime 
// whilst test prime <= ,/p: 
whi le( (( long) * tp) * { *tp ) <=p ) { 

i f ( ! ( p%*tp ) ) { // if current test prime divides (is factor of) p: 

P+=2; // try next odd number 

if(p<*pp) { // done when p overflows: 

*++pp=0 ; // terminate list 

return ; 

} 

tp=prime+l; _ // point to first odd test prime 

else ++tp ; // Else point to next test prime 

// no prime <= yp divides p so p must be prime: 

*++pp=p: // so store it next in the list 

P += 2; // try next odd number 

if(p<*pp) { // done when p overflows: 

*++pp=0 ; // terminate list 

return ; 

} 

1 

} 
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S (n ) . C 



void getSF ( u_long n, struct SF_struct *SF) 

{ 

u_int *pp; // ptr to current prime 

u_long r; // 'residue' of n remaining for factoring 

SF->sfk=0 ; // no. of distinct prime factors discovered 

r=n ; 

pp=prime; // point to start of prime list 



// whilst current prime <= ^r and prime list not exhausted: 
while! (( long) *pp ) * ( *pp ) <=r && *pp) { 

if(!(r%*pp)) { // if current prime is a factor of r: 

SF->sfp[SF->sfk]=*pp; // store current prime as next prime 

SF->sfa [SF->sf k] =1 ; // set its multiplicity to 1 

r/=*pp; // 'divide out' current prime 

whi le ( ! ( r% *pp ) ) { // while current prime factors r: 

SF->sfa [SF->sf k] ++ ; // increment multiplicity 

r/=*pp; // 'divide out' current prime 

} 

^ SF->sfk++; // increment count of distinct prime 

++PP ; // next prime 

} 



of SF 



factors 



if(r»>l) { // If n contains prime > 2~16: 

SF-> sf p [SF- > s f k] =r ; // store it as last prime of SF 

SF->sf a [SF->sf k] =1 ; // set its multiplicity to 1 

SF->sfk++; // increment count of distinct prime factors 
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S(n) .C 



u_long S(u_long n) 

{ 

struct SF_struct SF; // to store SF of n 

int sfi; // index of current term of SF of n 

u_long Sn ? // current guess at S(n) 

u_long x; // Stcurrent term of SF of n) where it might exceed 

// current value of Sn. 

if(n==l) return 0; // special case 

getSF(n, &SF); // calc, and store SF of n 

// First guess at S(n) is S(p~a), where p is the largest prime in the SF 
// of n and a is its multiplicity. This pre-empts the calculation of S(p~a) 
// for the remaining terms where, as is likely, p*a for these terms is <= 

// this initial guess (since S(p“a) <= p*a always): 
sf i=SF. sfk-1 ; 

Sn=Spa (SF.sfptsfi] ,SF.sfa[sfi] ) ; 

while(sfi>0) { // while more term(s): 

sfi--; // next term 

i f ( SF . sf p [sf i ] *SF . sf a Csf i ] >Sn ) { // if this term may have larger S(p~a): 

x=Spa ( SF . sfp [sf i ] , SF . sf a [sf i ] ) ; // calc, it 

if(x>Sn) Sn=x; // if new max., update Sn with it 

} , 

} 

return Sn; // That’s all folks! 



u_long Spa(u_long p, int a) 

{ 

// Refer to item 3) of the algorithm description in S(n).h. 

int c; // largest multiple of p such that f(p*c, p) <= a (eventually!) 

int z; // f(p*c, p) 

int m; // used to calc. no. of times p appears as factor of c 
if(a<=p) return p*a; 
c=a-2 ; 

z=f (p*c , p ) ; 
while(z>a) { 

// d in items 3.2. 3.1) and 3. 2. 3. 3) of algorithm description is implicit 
// here: 

2 — ; 
m=c-- ; 

whi 1 e ( ! (m%p ) ) [ // while p divides m: 

z — ; 

m/=p; // 'divide out' factor of p from m 

} 

} 

if(z<a) return p*(c+l); 
else return p*c; 
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S ( n ) - C 



int f ( int x, int p) 

{ 

int k=0? // count of appearance of prime p as a factor in the integers 

// from 1 to x. 

int xdp; // successive divisions of x by p 

xdp=x/p ; 
while (xdp>0) { 
k+=xdp ; 
xdp/=p ; 

} 

return k; 

} 
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for n«4294965296 to n - 4 29 4 9 6 7 Z 9 5 



P»«« 1 of 2 



SMARAMDACJ* FUNCTION, S( n ) , 
\S of n plus: 



4294965296 



4294965316 

4294965326 

4294965336 

4294965346 

4294965356 

4294965366 

4294965376 

4294965366 

4294965396 

4294965406 

4294965416 

4294965426 

4294965436 

4294965446 

4294965456 

4294965466 

4294965476 

4294965486 

4Z9496S496 

4294965506 

4294965516 

4294965526 

4294965536 

4294965546 

4294965556 

4294965586 

4294965576 

4294965566 

4294965596 

4294965606 

4294965616 

4294965626 

4294965636 

4294965646 

4294965656 

4294965666 

4294965676 

4294965686 

4294965696 ' 

4294965706 

4294965716 

4294965726 

4294965736 

4294965746 

4294965756 

4294965766 

4294965776 

4294965786 

4294965796 

4294965806 

4294965816 

4294965826 

4294965836 

4294965846 

4294965856 

4294965666 

4294965878 

4294965886 

4294965896 

4294965906 

4294965916 

4294965926 

4294965936 

4294965946 

4294965956 

4294965966 

4294965976 

4294965966 

4294965996 

4294966006 

4294966016 

4294966026 

4294966036 

4294966046 

4294966056 

4294966066 

4294966076 

4294966086 

4294966096 

4294966106 

42949661 1 6 

4294966126 

4294966136 

4294966146 

4294966156 

4294966166 

4294966176 

4294966 186 

4294966196 

4294966206 

4294966216 

4294966228 

4294966236 



l\ 8 


1 


2 


3 


i 268435331 


188887 


2739847 


60492469 


i 5531 


4294965307 


1194373 


101977 


i 1073741329 


4583741 


102261079 


27356467 


1 2147482663 


648101 


7255009 


2671 


i 21517 


11392481 


5520521 


2199163 


l 3919 


4294965347 


6911 


90407 


i 14969 


15859? 


3 109 


108907 


i 419 


30971 


4925419 


1431655123 


i 5527 


830587 


l 54573 


279857 


, 2147482693 


21367987 


3567247 


30034723 


357913783 


4040419 


35481 


27012361 


235289 


2347 


1 54807 


529393 


5179 


625997 


22138997 


4013 


816223 


§965117 


12109 


6040739 


65983 


1530091 


79536397 


2473 


178793 


§5579 


36097 


564607 


8329 


4294965457 


567067 


423191 


18275037 


14364433 


7615187 


1822217 


157 


1194041 


2147482739 


1935541 


102261083 


4234965487 


14639 


973253 


50951 


1 4387 


715827583 


390451409 


17903 


72031 


29483 


3804221 


3579 13793 


99862919 


4699087 


12236369 


2147482763 


252645031 


57193 


4294965529 


5835551 


26881 


88919 


3296213 


19979 


1 07077 


1073741387 


23333 


49481 17 


20355267 


1038937 


12569 


13174741 


1 431655189 


87108837 


397057 


16268809 


1618909 


3081037 


1431655193 


19701877 


139 


2843 


2549 


1073741399 


216941 


961703 


226050821 


154373 


118361 


536870701 


923053 


515231 


4294985817 


2957 


8961 1 


692513 


1828423 


48684409 


2459 


9177277 


682333 


2147482819 


1431655213 


2937733 


209623 


89476451 


119747 


20731 


26371 


147827 


4294965659 


9061109 


171 11417 


119557 


4603393 


1073741419 


28825273 


715827613 


4294965679 


195225713 


1431655223 


400949 


2038427 


* 2729 


58741 


1753 


477218411 


1059439 


594353 


1151 


75431 


220889 


1286303 


2147482859 


889411 


1716813 


849311 


601873 


8783161 


536870717 


423496573? 


179 


7517 


2147482873 


887573 


63161281 


650851 


2818219 


4294965757 


7877 


204522173 


1026031 


4294965767 


2311 


24793 


268435361 


1404961 


1162887 


603969 


4286393 


244157 


84407 


1811 


4337 


55823 


57649 


919 


23831 


41357 


24317 


116080157 


6761 


677113 


349241 


2165893 


40829 


2129 


357913819 


4799 


359231 


7883 


5355319 


4294965839 


1016247 


4294965847 


36599 


1431655263 


7129 


237173 


72211 


148102271 


107069 


1188127 


1073741467 


33353 


63113 


36749 


30678327? 


477218431 


2147482943 


4294365687 


2179 


2081903 


19289 


84121 


2147482949 


613566557 


6871 


364691 


1073741477 


110127331 


1206451 1 


75217 


167759 


16481 


258049 


3119 


1574401 


52181 


8947845? 


4294965937 


1621 


130150483 


7405113? 


263 


3933119 


4294965949 


5417 


38893387 


6221 


138547289 


1191061 


4294965967 


2591 


8699 


76695821 


4294965977 


.68389 


1321 


7430737 


367091 1 1 


26188817 


15259 


6907 


5956957 


33331 


4289 


36398017 


4294966007 


178958917 


60492479 


15101 


14737 


415133 


2579 


279511 


35495587 


258173 


254879 


5651271 1 


2235797 


65075243 


47197423 


236001 


163003 


134217689 


273617 


4519 


2213 


2147483023 


59069 


2147483033 


8 13566581 


32Z4449 


9695183 


103613 


477218453 


182129 


91382257 


79536409 


4294966087 


3701 


1431655363 


20543 


44278001 


45137 


4294966099 


2147483053 


683937 


1073741527 


35083 


131441 


3271109 : 


2147483059 


3329431 1 



2749 
536870767 
1993949 
4858559 
1 5675059 
481067 
126613 
1012009 
1063637 
14510021 
16007 
1121987 



73583 

521 

340519 



257863 

2147483069 

8873897 



4 

14316551 


5 

3226871 


6 

49941457 


7 

231547 


8 

337867 


9 

2798023 


5 1 74657 


66569 


987 


4294965313 


1 13025403 


286331021 


1602599 


4919777 


2147462661 


22967729 


211 7833 


24542659 


57427 


4294965331 


3741253 


75577 


85889 


858993067 


12632251 


40139863 


2286969 


535613 


979691 


462 1 


6364 1 


14316551 17 


1228537 


2545919 


715827559 


1 493 


2339 


4294965361 


2 1 4748268 1 


23719 


118423 


38629 


195137 


1189741 l 


35731378 l 


104755253 


64 1 


1307 


1 199711 


204522161 


39107 


4294965383 


6628033 


858993077 


7535027 


35449 


268435337 


5749619 


2147482697 


86357 


1352257 


377513 


70867 


2634948 1 


51 1 549 


977239 


5167 


14461 163 


672349 


4294965413 


34301 


24683 


23860919 


15107 


17747791 


204522163 


39727 


10105801 


321239 


28793 


5772803 


746561 


5003 


1453457 


13421787 


1551089 


308857 


2617 


28871 


858993089 


26633103 


11731 


121313 


1399? 


2465537 


2207 


70019 


4294965461 


415937 


138547273 


536870683 


168331 


40699 


88873 


163681 


4001 


5549051 


6079 


18233 


2203 1 


7 18943 


84869 


1073741371 


8863 


18673763 


2754349 


7253 


30509 


2866 1 


137593 


1227133 


635759 


16909313 


837519 


294337 


72679 


255197 


Z 0 1 7 


103583 


7033 


306783251 


188417 


53687069 


832519 


71582758? 


142799 


1 07374138 1 


1347297 


364509 


103391 


32789 


6381821 


12132671 


122713301 


2468371 


6529 


230664 1 


1 1831861 


64319 


106747 


2321603 


822549 


89478449 


179623 


17599 


891997 


28057 


2293 


756 1 


77933 


119304599 


78090283 


20452217 


19709 


51829 


20507 


298801 


103 1 


11302541 


4294965561 


10374013 


12781 


2949839 


337 


429496559 


9739151 


8801159 


72796027 


2970239 


50529007 


198641 


4294965601 


2147462801 


33294307 


40459 


533536 1 



9283 


35527 


307 


4294965613 


78487 


21859 


5147 


7435577 


214748281 1 


330381971 


4409 


3061 


251809 


9137 


170327 


1302669 


2147482817 


3 l 93283 


2551 


4294965641 


784039 


702? 


4549 


40153 


85899313 


62245879 


100547 


6400843 


160033 


40961 


1745921 


16582879 


118757 


28071671 


20929 


858993133 


429496567 


4294965671 


59652301 


4294965673 


93368819 


3011 


101873 


18013 


505409 


4294965683 


1360889 


12347 


47721841 


4294965691 


130261 


45613 


1 130254 13 


7890671 


42949657 


3373 


94399 


1305460? 


201907 


286331047 


613 568 53 


13287 


72727 


2238127 


1831 


858993143 


1556147 


4294965721 


118273 


249287 


153391633 


29153 


70099 


494413 


23057 


4294965733 


2147482867 


859853 


16519039 


130150477 


1689601 


7309 


6701 


1553333 


121843 


226050823 


29173 


142837 


2147482877 


6270023 


1877721 


18755309 


89669 


82183 


1 07374144 1 


4457 


7229 


477218419 


1000691 


19792469 


65075239 


65497 


4933 


12377423 


836573 


42223 


536870723 


16943 


3539 


41638697 


7433 


4294965793 


306783271 


286331053 


21474829 


7052489 


1 02559 


S 5 6 5 5 9 


1451 


8753 


143165527 


3S008S47 


1 0627 


1431655271 


17471 


122713309 


57527 


4294965621 


20789 


81037091 


20719 


19088737 


38669 


651641 


536870729 


1335629 


238609213 


4999 


1626881 


4294965841 


2147462921 


4866 1 


1757351 


99409 


721843 


954649 


32743 


346 1 17 


4303573 


733 


214748293 


11096103 


309391 


4864061 


25565273 


7933 


407879 


613566553 


268435367 


65173 


2147482937 


1 5959 


107374147 


186737647 


28031 


70309 


1091 


5843491 


2161 


366997 


6869 


620929 


6949783 


78090289 


14316553 


252645053 


2147482951 


9833 


1009 


858993181 


429496591 


429496591 1 


87767 


7829 


7588279 


944987 


1 03643 


405683 


25799 


125941 


351931 


171798637 


181913 


97553 


831713 


2243 


43826183 


3886847 


214748297 


49003 


839703 


1290167 


536870743 


19927 


4651 


79259 


3947579 


25718359 


1494421 


122713313 


627919 


30034727 


306783283 


1431655321 


241453 


858993193 


991909 


4294965971 


13597 


311749 


1390253 


457 


690509 


1431655327 


911 


74471 


91867 


80513 


229 


61333 


1 17709 


528286 1 


5276371 


858993199 


131 


4294966001 


71 5827667 


1322749 


38053 


32497 


13159 


243769 


107374 1 503 


471301 


1318283 


1340083 


5939 


40031 


1 9701673 


6537239 


529981 


192 19 


2225371 


116080163 


1439 


3637 


358691 


11003 


110581 


320927 


93368827 


4294966043 


357913837 


1 1 2S769 


9544369 


3256229 


3343 


238649 


126322531 


277363 


2788339 


54366659 


24683713 


351263 


28627 


16 231 


10475527 


3853 


940223 


4294966073 


11987 


171798643 


1033 


12521767 


303617 


46182431 


71 1 0871 


6879 1 


4013987 


11392483 


357913841 


1 081563 


8803 


13634813 


829 


85733 


306783293 


340573 


2657 


1433 


233549 


4554577 


7064069 


167269 


69273647 


78090293 


22171 


4Z949661Z1 


1403583? 


2633333 


1823 


1 1 453243 



4294966129 
252645067 
9697 



674249 

1213 

40009 



1759 

1328477 

8641783 



242873 
58040083 
4 1 S3 



665783 
792283 
4294966153 



1 122943 
4793489 
2147483077 



2340581 

27709459 

286331077 



4294966246 2147483123 



4294988256 

4294966266 

4294986276 

4294966286 

4294966296 

4294966306 

4294966316 

4294966326 

4294966336 

4294968346 

4294966356 

4294966366 

4294966376 

4294966386 

4294966396 

4294966406 



36347913 
163543 
85968 1 
16519101 1 
304867 
251 197 
26188819 
1571 
17159 
195225743 
56693 
263571 
26249 
3851 
32783 
1 16099 



186737659 


77647 


390451469 


53687077 


5168431 


288679 


4294966163 


578213 


122713319 


1 55801 


905347 


251 


13015049 


82139 


16087 


666 1 


9629969 


171798647 


4294988 1 77 


74051141 


2237 


9336883 


390451471 


55063669 


817933 


97987 


86531 


4294986187 


18837571 


80737 


119737 


21367991 


24403217 


14385039 


191 


365063 


47279 


2008871 


7548271 


1022611 


352133 


6317 


183 13 


2562629 


19976587 


2113 


516ZZ19 


57847 


1633067 


2245147 


525571 


6719 


21503 


286331081 


4294988217 


43481 


237409 


10853 


1748053 


30246241 


1361 


5283439 


15618059 


890887 


9787 


446323 


193673 


4294966231 


536870779 


1739557 


4029049 


512671 


4294966237 


95287 


1130071 


246271 


36239 


23743 


4294966243 


1 073741 56 1 


15070057 


566693 


482421 


28676609 


3307 


1431655417 


1 109 


1 16080169 


763957 


9323 


477218473 


171401 


12744707 


13463 


5501 


16393001 


29671 


SZ 1 


29620457 


355573 


1073741567 


361673 


61356661 


96263 


683 


2217329 


2147483137 


1 79 


87652373 


715827713 


2363 


9761287 


4073 


6691 


35537 


58 169 


858993257 


15391 


1 43779 


16417 


2663 


311 


63161269 


6721387 


74051 143 


37347533 


4294966297 


3883333 


332711 


26141 


72736033 


380557 


35747 


1254371 


4090444 1 


2458481 


357913859 


127663 


919693 


477218479 


2753 


390451483 


15230377 


31547 


5940479 


108191 


6074917 


5051 


226050859 


139801 


49367429 


160021 1 


197243 


1 16080171 


538870731 


159072827 


430789 


t 94 8 9 


174337 


613566619 


49941469 


397 


4294966337 


1411889 


325697 


1614649 


2706343 


2147483171 


252645079 


178956931 


858993269 


68174069 


77813 


39403361 


28633109 


135347 


5209 


54 133 


200381 


856993271 


876 1 


2147483179 


75350287 


107053 


146681 


23091217 


1 762393 


1367 


25733 


429496636? 


58 1 029 


401887 


429496637 


1431655457 


117877 


4294966373 


170557 


70 1219 


1907 


1 13025431 


390451489 


14251 


908603 


108827 


159072829 


268435399 


73727 


186737669 


151637 


2464123 


273543 


33818633 


890333 


1057613 


126322541 


1999 1 


1013681 


17459213 


399941 


1342177 


5557 


48761 


109001 


3579 13867 


858993251 


1431655469 


6047 


41698703 


3407 


22067 


9761 2873 


149333 


1679033 


655993283 
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SMARAJ*OAO« FUNCTION. S< n > , for n-4294965296 to n-4Z94967Z95 
\S of n plus: 

ft \ 0 123458759 
4294986418 1S77 813566631 48677 8066449 251 30899039 4604213 390451493 7213 57266219 
4294966426 1564081 4294986427 7017919 95569 2939 204522211 577 389649 12341 S36763 
4294968436 37571 17674739 809 15913 Z81251 4294966441 614093 1431655481 41539 656221 
4294966446 3911827 4294966447 1172207 1137137 965181 2669339 78301 3929521 199 5039 
4294966436 28861 288427 6567227 172787 202021 19611719 126322543 99882941 2485513 858993293 
4294966466 306783319 130150499 2848121 17970571 143165549 12763 6468323 1949 2147483237 13381 
4294986476 38891 4294966477 3308911 7001 881 330382037 321143 2176871 1073741621 2951867 
4294966486 89273653 613566641 178956937 3455323 429496649 159072833 56512717 6806603 7866239 858993299 
4294966496 7895139 49367431 2147483249 4999961 86787 813568643 5413 1431655501 536870813 6151 
4294966506 234467 416623 1153 477218501 81083 1868189 443 390487 1256573 26407 
4294966516 306871 49891 1951 124193 8259351 345727 3984199 9565627 10631 7469507 
4294966526 2251 477218503 16777213 613566647 443237 1210873 4493 32611 29417579 40123 
4294966536 8521759 482093 2147483269 160S60I7 396347 2307881 715827757 613566649 149713 56509 
4294966546 165191021 252545031 16193 3023 12271333 1431855517 538870819 4294966553 238609253 32443 
4294986558 73883 1129957 8910719 330382043 84373 119083 52377841 617359 530787 16207421 
4294966588 34781 85563 483219 511123 201547 25577 5107 186163 846199 364753 
4294966576 1405421 293 37633 138547309 214748329 21577 22139003 4294966583 6170929 9439487 
4294966586 58040089 75350291 133169 12739 47721851 4294966591 222953 21367993 298303 179593 
4294966596 383 4441 9717119 888627 32587 465781 20369 Z08261 1073741851 286331107 
4294966608 15661 58835159 29828157 22253713 1523 110127349 5088823 7993 9061111 490573 
4294966616 870131 105323 764501 4294966619 6551 1693 21262211 92347 11719 16901 
4294966626 79536419 45751 12064513 23469763 1297573 312839 719 8597 306783331 2219621 
4294966636 873871 330382049 7639 4294966639 5621 15671 74051149 14549 153151 4567 
4294966646 2147483323 1281697 1080223 284831 2202547 4294966651 1073741683 799361 195225757 15877 
4294966656 17449 4294966657 4121849 139361 86627 4Z94966661 20063 756023 63997 26030101 
4294966666 556467 4294966667 152239 214309 429496667 1237 8659207 1033189 715827779 120223 
4294966676 757 8251771 49941473 18979 8053 24626397 86477 12030719 1073741871 5867 
4294966686 17959 3*0451517 2383 36709117 21401 109451 357913891 54366667 1028981 286331113 
4294966696 536870837 16582883 21691751 109961 87119 75350293 38569 138547313 12782639 218629 
4294968706 2147483353 975907 4492643 16976153 143165557 2280917 536870639 1431655571 701 2111 
4294966716 119304631 3607 919 20164163 431 2194 462689 69371 20259277 4241 
4294986726 52377643 64103981 13765919 208689 144563 64349 397 889043 4079 924643 
4294966736 1020887 30460757 11864549 32292983 71582779 40129 195225761 559459 36209 658993349 
4294966746 115811 18041 38653 1431655583 17179867 4141723 14913079 69463 79841 35671 
4294966756 1073741689 5811863 1299143 198391 4789 89417 46411 15176561 756689 27709463 
4294966766 321047 191 27751 4294966769 47721853 118080183 107021 24265349 306783341 15618061 
4294966776 1844917 9133 99223 477218531 16519103 2387419 715827797 6938557 178481 154523 
4294966788 1490273 148102303 34019 85199 7459 562979 536870849 330382061 10663997 78311 
4294966796 1789 43383503 2147483399 17821439 3579139 93229 93503 1431655601 4481 317323 
4294966806 761 99662949 32173 204522229 423496681 1072937 32779 4294966813 4021 88129 
4294966816. 3257 1516049 359171 1766749 214746341 13899569 26561 2907899 25679 61291 
4294966826 29417581 13229 9502139 4294966829 35081 91382273 349981 477218537 27059 858993367 
4294966636 48623 613566691 15449521 88609 349753 38351 20347 19259941 275369 1579 

liuima iM 7 i??si 7 imi *miu ..tm 7 !iai ??s? s ji?? 

4294966686 191449 2999 1073741717 2909 293 330382067 1345541 1327 6H993 2290649 
4294966676 1073741719 4294966677 5076793 911689 26843543 1431655627 58040093 138053 27531839 221219 
4294966886 3067S3349 338213 23633 1931 3491843 8311 1073741723 911 50647 9651611 
4294966896 1489 330382069 1246363 154423 6135667 113783 715827817 59791 22193 95443709 
42949669«€ 5639 178207 357913909 4294966909 2029 1431655637 70051 13687 2797 858993383 
4294966918 7307 10928609 82339 641327 1154561 36092159 6555711 453199 3271 36191 
4294966926 715827621 4294966927 572357 8263 5436667 104755291 57943 12744709 35204647 21839 
4294966936 16121 2403451 42107519 186737693 19522577 8297 306783353 4294966943 9151 858993389 
4294966946 20069939 5990191 56512723 47197439 53923 9817 536670869 1431655651 2147483477 683 
4294966956 1381907 353699 23067 397351 36497 447439 217643 87652367 1073741741 393853 
4294966966 70067 226050893 59852319 97169 7621 1431655657 449 390451543 715827829 12113 
4294966976 67108859 163487 2147483489 2S0791 347 4294966981 20297 1431655661 78341 41177 
4294966966 12558383 569851 602887 1588963 98081 3427747 323027 18341 2147483497 193 
4294966998 119159 4294966997 26107 16330673 4294967 110127359 4106087 554977 13256071 6458597 
4294967006 11483869 65519 5835553 1597 143165567 22453 21913097 375467 12707003 $58993403 
4294967016 38149 795217 14412641 150053 180007 365747 5549053 17401 14128181 10141 
4294967026 23917 18541 32537629 4294987029 11608019 5683 536870879 3764213 3392549 858993407 
4294967036 2445881 1413283 8355967 20549 241 10710641 96703 9851 34636831 11767033 
4294967046 18111 39163 141319 477218561 5009 21582749 357913921 37201 2879 1228889 
4294967056 5347 252645121 83459 1755197 214748353 7121 4007 768193 564533 19976591 
4294967066 15761 132967 2599859 2408843 143165569 65551 12201611 16086019 1779191 791699 
4294967076 119304641 23761 5250571 110127361 107374177 417839 102281121 390451553 342283 12917 
4294987086 2147463543 4294967087 1132639 613566727 429496709 21467 9209 11332367 21691753 395303 
4294967096 76695841 1431655699 2147463549 96553 362 181967 66977 68174081 5197 9199 
4294967106 1583 2135737 2609 1431655703 1252177 4294967111 182423 88463 675097 286331141 
4294967116 63841 10058471 754297 15053 236507 102871 12859183 17321 81031 14753 
4294967126 2147483563 21841 $5639 161471 463319 12899 1073741783 60021 186689 129347 
4Z94967136 13147 138547327 1467863 159072857 214748357 1163 2476913 4294967143 7321 576119 
4294967146 93368851 13990121 2251031 2844249 85899343 77249 16381 33613639 70853 515293 
4294967156 1073741789 1932059 2147483579 5156023 647 4294967161 34919 4886197 381707 858993433 
4294967166 1262483 3229 1801 179383 184571 651839 357913931 9739 2147483587 7283 
4294967176 874*1 104755297 317159 62929 34981 3583 195225781 14787 6197 44171 
4294967186 152293 204522247 5449 4294967189 143165573 82791 23342213 29333 7134497 45210181 
4294967196 2521 4294967197 1296007 46182443 31033 11576731 23S609Z89 26029 97612891 286331147 
4234967206 1809169 91382281 18539 488009 25264513 477218579 292493 75167 85981 13109 
4294967216 3780781 1431655739 52377649 11959 25411 39403369 813749 3089 1726273 231223 
4294967226 22037 252645131 10631107 69779 13854733 4294967231 8191 3557 318949 22025473 
4294967236 417961 825479 Z9389 70409299 23327 1431655747 145661 619139 45083 69001 
4294967246 26111 65543 1877171 15966421 4111 1664071 1527371 1431655751 *62097 4259 
4294967256 59652323 *2613 2147483629 4817 32771 281 19867 64103989 262657 95443717 
4294967266 367531 37549 29587 35609 115861 3449 536870909 395959 6116167 171798691 
4294967276 361651 1861711 18046081 4294967279 127 8951 795659 1226783 46684427 122713351 
4294967286 715827681 65539 2089 1039 22605091 4294967291 331 464773 2147483647 65537 




Timing of S(n).c Module for Calculation of Smarandache Function, version 1 



Time taken to calculate S(n) depends on how easy it is to factor n. Less 
time is required if n has "small" prime factors. So, in the following table, 
the values of n shown are the mid-points of ranges (n-99 thru n+99). Times 
shown are for calculating S(n) for all integers in each range 10 times over: 



n 


time (ms) 


100 


268 


200 


308 


400 


345 


800 


387 


1600 


432 


3200 


490 


6400 


571 


12800 


661 


25600 


766 


51200 


919 


102400 


2450 


204800 


4036 


409600 


5670 


819200 


7977 


1638400 


10423 


3276800 


13004 


6553600 


16302 


13107200 


23438 


26214400 


29642 


52428800 


37011 


104857600 


50330 


209715200 


62363 


419430400 


77888 


838860800 


108179 


1677721600 


158480 


3355443200 


159111 
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Timing of S(n).c Module for Calculation of Smarandache Function, version 2 



"Time to n" is the 
"Time add. " is the 
<= current n. All 



time taken to calculate S(n) for all n <= that shown, 
time taken to calculate S(n) for all n > previous n and 
times are in milliseconds (as per version 1): 



n 



50000 
100000 
150000 
200000 
250000 
300000 
350000 
400000 
450000 
500000 
550000 
600000 
650000 
7 0 00 0 0 
750000 
800000 
850000 
9Q0000 
950000 
1000000 
1050000 
1100000 
1150000 
1200000 
1250000 
1300000 
1350000 
1400000 
1450000 
1500000 
1550000 
1600000 
1650000 
1700000 
1750000 
1800000 
1850000 
1900000 
1950000 
2000000 
2050000 
2100000 
2150000 
2200000 
2250000 
2300000 
2350000 
2400000 
2450000 
2500000 



Time to n 



18223 

66763 

139191 

229634 

335252 

452539 

579419 

715146 

859963 

1012335 

1171221 

1336899 

1508825 

1686808 

1870023 

2058983 

2252457 

2450892 

2653620 

2860734 

3072049 

3288502 

3509106 

3733965 

3962171 

4194158 

4429331 

4668560 

4910513 

5155601 

5404652 

5656512 

5911306 

6169686 

6431383 

6696172 

6963206 

7232974 

7505412 

7779763 

8056579 

8336442 

8620053 

8905641 

9194727 

9486449 

9780105 

10076920 

10375202 

10676383 



Time add. 



18223 

48540 

72428 

90443 

105618 

117287 

126880 

135727 

144816 

152372 

158886 

165678 

171927 

177983 

183215 

188961 

193473 

198435 

202728 

207115 

211314 

216454 

220603 

224860 

228206 

231987 

235173 

239229 

241953 

245088 

249051 

251859 

254794 

258380 

261697 

264789 

267034 

269768 

272438 

274351 

276816 

279863 

283611 

285588 

289086 

291722 

293655 

296815 

298282 

301181 
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